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ABSTRACT 



A sudden increase in stellar luminosity may lead to the ejection of a large part of any optically thin gas orbiting the star. 
Test particles in circular orbits will become unbound, and will escape to infinity (if radiation drag is neglected), when 
the luminosity changes from zero to at least one half the Eddington value, or more generally, from L to (I/Edd + L) /2 
or more. Conversely, a decrease in luminosity will lead to the tightening of orbits of optically thin fluid. Even a modest 
fluctuation of luminosity of accreting neutron stars or black holes is expected to lead to substantial coronal heating. 
Luminosity fluctuations may thus account for the high temperatures of the X-ray corona in accreting black hole and 
neutron star systems. 
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1. Introduction 



Low mass X-ray binaries (LMXBs) are bright X-ray sources 
whose luminosity is thought to be powered by accretion 
from a binary companion onto a neutron star or a black 
hole. Typically, each source can be found in one of two 
or more distinct spectral states, some of which include a 
hard X-ray power law component intcrpctcd as radiation 
from a hot (~ 10 2 kcV) corona. Each state has an associ- 
ated characteristic variability on many timescales. Nearly 
all sources exhibit quasi-random excursions of luminos- 
ity. In this Letter I examine some of the consequences of 
rapid variations of luminosity, particularly for optically thin 
Keplcrian flow, which I identify with the X-ray corona. 



Since the paper of Walker & Meszaros (1989) 



it is un- 
derstood that high luminosity of an accreting neutron star 
will be associated with removal of angular momentum from 
the optically thin flow via the Poynting-Robertson drag. In 
fact, it seems to be a part of the lore of neutron star as- 
trophysics that the primary effect of high luminosity is to 
bring matter down from its nearly-circular orbits towards 
the inner parts of the accretion disk, and perhaps even onto 
the stellar surface. This is ideed true for optically thin flow 
in a steady radiation field. However, a rapid change in lu- 
minosity may have the opposite effect and, as I show below, 
under certain conditions may lead to ejection of the opti- 
cally thin matter. 



An even more important effect will occur when the lu- 
minosity undergoes a small but rapid change. The optically 
thin fluid will then continue its orbital motion in very nearly 
the same circular orbits, but will undergo strong heating, 
caused by dissipation of excess kinetic energy. This may 
explain the presence of hot X-ray coronae in some spectral 
states of LMXBs. 



2. Keplerian orbits 

In Newtonian physics, computing the response of an or- 
biting test particle to a sudden change in luminosity is a 
simple exercise in orbital mechanics. Consider orbits in a 
spherically symmetric potential 



V(r) 



GM 



(1) 



As shown by Newton these are conic sections with eccen- 
tricity 



e = \ 1 



2EP 
G 2 AP 



(2) 



where E and I are the specific (per unit mass) energy and 
angular momentum of the orbiting particle. The specific ki- 
netic energy and angular momentum of a particle in circular 
orbit of radius ro are 



GM 



h = x/GM 



r , 



and the total specific energy is 



E 



GM 
~2nT' 



(3) 



(4) 



Suppose that the energy of the particle is changed to E, 
with no change in its angular momentum. The particle will 
now move in an orbit with eccentricity 



= y/1 + E/ Kq 



(5) 



In particular, E ~ Eq =>■ e = 0, i.e., the original orbital en- 
ergy corresponds to a circular orbit; while E = =>■ e = 1, 
i.e., a marginally unbound particle moves in a parabolic or- 
bit. It is well known that E > corresponds to an unbound 
hyperbolic orbit. As we will see, the effect of an impulsive 
change of luminosity on test particles in circular orbits is 
equivalent to a change of the particle energy. 
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3. Luminosity effects 

Optically thin hydrogen plasma suffers radiation pressure 
forces proportional to the Thomson cross-section and the 
radiative flux. At Eddington luminosity, L^d, the radiative 
force balances gravity exactly at any radius — in Newtonian 
physics of a spherically symmetric source both the radiative 
flux L/(4:Trr 2 ), and the force of gravity —GM/r 2 drop off 
inversely with the square of the radial distance. Hence, at 
lower luminosity the presence of radiation pressure is equiv- 
alent to a proportional reduction of the gravitational mass 
by irradiation: GM may simply be replaced by GM(1 — A), 
with A = L/L^dd, so that eqs. (Q]), ([2} become 



V(r) = - 



GM(\ - A) 
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(7) 



The specific kinetic energy in circular orbit at ro for a source 
with luminosity L and true gravitational mass M is 



GM{1 - A) 
2r ' 



(8) 



Upon an impulsive change of stellar luminosity from 
zero to L, a particle that has been travelling travelling in 
circular orbit at ro would conserve its kinetic energy and 
angular momentum, K , Iq of Eq. ([3]), so that its new orbit 
would be described by the following specific energy and 
angular momentum: 



E = 



GM GM(1 - A) 
2r r 



l = l = y/GM 



ro, 



(9) 



and because of the excess energy, K — K\ = XGM/(2r ), 
this orbit would no longer be circular. The eccentricity of 
the orbit is, from Eq. (0, 



A 



1-A 



(10) 



In particular, A > 1/2 e > 1, i.e., if the luminosity is 
now one half of the Eddington value, or greater, the par- 
ticle is in an unbound orbit. If, on the other hand, the in- 
crease was from zero to less than one half of the Eddington 
luminosity, the particle remains in a bound elliptic orbit, 
A<l/2=>e<l, with an increased semi- major axis 
a = r (l - A)/(l - 2A). 



4. Coronal ejection 

In reality, the fluctuations of luminosity in LMXBs do not 
occur between L = and L / 0, but rather between some 
initial L\ ^ and some final value L ^ 0. Generalizing the 
derivation from the previous section, we start with circular 
orbits at luminosity L\, with Ai = Li/LEdd, the specific 
kinetic energy and angular momentum in circular orbits in 
the optically thin region being given by 



GM (1 - Ai) 



h = VGM(l-Xi)r. (11) 



Following an impulsive change of luminosity to L, the new 
orbits arc defined by a new potential and the same values 
of kinetic energy and angular momentum: 

V{r) = - GM{1 ~ X \ K = K U l = h , (12) 



where, as before, A = L/ Z/Edd- The eccentricity of the new 
orbits is 

„ |A-Ax| 



1 - A 



(13) 



Identifying the corona with the optically thin region, 
the condition for its ejection (e > 1) now becomes: 



i.e., 



1- A< (1- Ai)/2 



^Edd — L < (Lsdd — L\)/2, 



(L Edd + ii)/2 < L. 



(14) 



(15) 



(16) 



Thus, the condition for coronal ejection is that the luminos- 
ity reduces its distance to the Eddington value by at least 
a factor of two. Clearly, the closer the initial value L\ is to 
^Edd, the lower the fractional value of luminosity increase 
necessary for the ejection of the corona. E.g., an increase 
from 0.8iEdd to 0.9LEdd corresponds to a fluctuation of less 
than 13%, and this is sufficient to clear out the corona. 



5. Coronal heating 

Now we turn to changes of luminosity that do not lead to 
coronal ejection. This could be because the luminosity has 
increased (from L\) by less than (L^dd — L\)/2, thus violat- 
ing condition (|16p , or because the luminosity has decreased. 
In any case, the new test particle orbit is an ellipse with 
semi-major axis ro(l — A)/|l + Ai — 2A|, in the notation 
of the previous Sections. The periastron is at radius ro for 
A > Ai , and at 



= r (l - A x )/(1 + Ai - 2A), 



(17) 



for A < Ai. For a small decrease of luminosity, the latter 
value becomes r_ w r [l — 2(Ai — A)/(l — Ai)]. 

Naturally, the state of the corona cannot be described 
by test particle orbits, as those of different particles would 
intersect. It is clear that some energy will be dissipated, 
and the fluid orbits will circularize. The ultimate outcome 
is difficult to predict, because redistribution of angular mo- 
mentum may occur. However, it seems clear that the dis- 
sipated energy will be on the order of (Ai — A) 2 times the 
virial energy of the corona. 

Consider the evolution of a geometrically thin annulus, 
initially orbiting at ro, under a small impulsive change of 
luminosity of the central star, \SX\ = |A — Ai| << 1, and 
subsequent dissipation of energy. If the annulus conserves 
its angular momentum, it will settle down in a circular orbit 
at n = ro(l — Ai)/(1 — A) after dissipating an amount of 
energy equal to 

(SX) 2 GM 



1 - Ai 2r 



(18) 
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6. Discussion 

The observed accreting neutron stars and black holes are 
usually quite luminous and are typically variable in time. 
It seems important to explore the consequences of rapid 
changes of luminosity, whose magnitude may be a large 
fraction of the Eddington value. 

The relevance of Poynting- Robertson drag to neutron 



stars was first pointed out by Walker & Meszaros (1989) 
While there is no doubt that high luminosity is accom 
panied by radiation drag, which at least in the optically 
thin regime in a steady source will eventually remove ki- 
netic energy and angular momentum of orbiting matter, 
this process takes time. If the luminosity undergoes a rapid 
change, the first and immediate response of orbiting matter 
is to change its trajectory. 

By considering Newtonian orbital mechanics, I have 
shown that an impulsive increase of central luminosity 
of sufficiently high magnitude, L — L\ > (£Edd — Li)/2, 
may lead to an ejection of the optically thin corona on 
a dynamical timescale. This may have an application to 
X-ray spectral-state changes of black holes and neutron 
stars, and may be of some importance in X-ray bursts. 
Inclusion of radiation drag requires numerical computa- 
tions, and is postponed till another paper, where a fully 
general relativistic discussion of the problem will be pre- 



sented |Stahl et al. (2012) 

The coronal response to small changes in luminosity was 
considered in Section [SJ It seems inevitable that relatively 
minor changes in the luminosity of the central source lead 
to substantial energy dissipation. This could be the as yet 
unexplained mechanism of coronal heating. The estimated 
magnitude of the effect is rather large. The dissipated en- 
ergy at 50 Schwarzschild radii caused by a single excursion 
in luminosity of 10% Eddington (O.lLEdd) is on the order 
of 10 -4 of the coronal rest mass, i.e., it corresponds to a 
temperature of ~ 100 keV. At 5 Schwarzschild radii the 
same result will be obtained by a fluctuation of only 3% 
Eddington luminosity. 

In passing, we note that for a steady source of luminosity 
A in Eddington units, the orbital frequency in the optically 
thin regime is modified by a factor of y/1 — A. Care must 
be taken to account for this effect when interpreting the 
rcdshift/blucshift of any spectral features in terms of orbital 
motion. 



The results of Sections 2-4 were presented on October 

28, 2011 at the LOFT Science Meeting in Amsterdam. 

http://www.isdc.unige.ch/loft/index.php/meetings/loft-science-meeting 
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